AD-A213  137 


Naval  Research  Laboratory 


Washington  DC  20375-5000 


NRL  Memorandum  Report  6549 


DTJC  Fill  cr-*-: 

A  Solvable  Self-Similar  Model  of  the 
Sausage  Instability  in  a  Resistive  Z-Pinch 


DTiC 

k  OCT  0  51989  4 


Martin  Lam  pe 

Beam  Physics  Branch 
Plasma  Physics  Division 


September  20,  1 989 


Approved  for  public  release;  distribution  unlimited 


89  1  0  4  05  4 


SECURITY  CLASS-FCat  ON  o 


REPORT  DOCUMENTATION  PAGE 


a  restrict. ve  markings 


Form  Approved 
OMB  No  0704-0188 


la  REPORT  SECURITY  C. ASS. F  CAT  ON 

UNCLASSIFIED 


a  SECURITY  CLASSiF, CAT  ON  AuTh OR 


2b  DECLASSIFICATION!  DO  WAG  RAD  NG  SC-EDo-E 


4  PERFORMING  ORGANIZATION  REPORT  NUM3ER;Si 

NRL  Memorandum  Report  6549 


6a  NAME  OF  PERFORMING  ORGANiZAT .ON  60  OFFICE  SVV80L  7a  NAME  OF  MONITORING  ORGANIZATION 

(If  applicable) 


3  DISTRIBUTION  •  AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
uni imi ted . 


5  MONITORING  ORGANIZATION  REPORT  NUM3E R(Si 


Naval  Research  Laboratory 


6c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Washington,  DC  20375-5000 


Code  4790 


7b  ADDRESS  (C/fy,  State  and  ZIP  Code) 


8b  OFFICE  SYMBOL  9  PROCUREMENT  NSTRuMENT  DENTIFiCAT ,ON  NUMBER 

(If  applicable) 


'0  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 
Arlington,  VA  22217 


1 1  Title  (Include  Security  Classification) 

A  Solvable  Self-Similar  Model  of  the  Sausage  Instability  in  a  Resistive  Z-Pinch 


PROJECT 

r  ASK 

NO 

NO 

ONR 

WORK  UNiT 
ACCESS.ON  NO 


12  PERSONAL  AUTHOR(S) 

Lanroe,  Martin 


13a  TYPE  OF  REPORT 

Interim 


16  supplementary  notation 


13b  TIME  COVERED 
FROM  _  tq 


14  DATE  OF  REPORT  (Year.  Month.  Day)  15  RAGE  COUN' 

1989  September  20  38 


COSATI  COOES 


18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  numoer) 


SUB-GROUP 


Z-pinch 

Pinch 


Sausage  instability1' 
MHD  instability 


19  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

A  solvable  model  is  developed  for  the  linearized  sausage  mode  within 
the  context  of  resistive  MHD.  The  model  is  based  on  the  assumption  that 
the  fluid  motion  of  the  plasma  is  self-similar ,  as  well  as  several 
assumptions  pertinent  to  the  long-wavelength  limit.  The  perturbations  to 
the  magnetic  field  are  not  assumed  to  be  self-similar,  but  rather  are 
calculated.  Effects  arising  from  time  dependences  of  the  equilibrium, 
e.g.,  current  rising  as  t  ,  ohmic  heating,  and  time  variation  of  the  pinch 
radius,  are  included  in  the  analysis.  The  formalism  appears  to  provide  a 
good  representation  of  those  modes  that  involve  coherent  sausage  distortion 
of  the  entire  cross  section  of  the  pinch,  but  excludes  modes  that  are 
localized  radially,  and  higher  radial  eigenmodes.  For  this  and  other 
reasons,  it  is  expected  that  the  model  underestimates  the  maximum 

(Continues) 


20  DISTRIBUTION 7 AVAILABILITY  OF  ABSTRACT  |21  ABSTRACT  SECURITY  CLASSIFICATION 

QUNCLASSIFIED/UNUMITED  □  SAME  AS  RPT  □  OTIC  USERS  I  UN 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

Martin  Lanroe  (2021  767-4041  Code  4?9Z 


DD  Form  1473,  JUN  86  Previous  editions  are  obsolete  _ security  classification  qf  'mis  £age 

S/M  0102-LF -01  'i -6603 


i 


5iF:CA7'C\  OF 


19.  ABSTRACTS  (Continued) 

instability  growth  rates,  but  is  reasonable  for  global  sausage  modes.  The 
net  effect  of  resistivity  and  time  variation  of  the  equilibrium  is  to 
decrease  the  growth  rate  if  a  <  1,  but  never  by  more  than  a  factor  of  about 
two.  The  effect  is  to  increase  the  growth  rate  if  a  >  1. 


00  Form  1473,  JUN  86  (Reverse) 


CONTENTS 

I.  INTRODUCTION  .  1 

H.  SELF-SIMILAR  EQUILIBRIUM  .  4 

m.  LINEARIZED  PERTURBATION  ANALYSIS  .  8 

IV.  SOLUTION  OF  THE  PERTURBATION  EQUATIONS  .  16 

CONCLUSIONS  . '  22 

ACKNOWLEDGMENTS  .  23 

REFERENCES  .  24 

DISTRIBUTION  LIST  .  29 


J 


A  SOLVABLE  SELF-SIMILAR  MODEL  OF  THE 
SAUSAGE  INSTABILITY  IN  A  RESISTIVE  Z-PINCH 


I .  Introduction 

Since  the  early  days  of  controlled  fusion  research,  it  has  been 

generally  accepted  that  simple  Z-pinches  are  violently  unstable  to  both  the 

sausage  and  kink  modes.  Specifically,  ideal  MHD  unequivocally  predicts 

instability  for  a  time-independent  pinch  surrounded  by  vacuum,  with  growth 

rate  of  the  order  of  kv.  for  ka  <1,  where  k  is  the  axial  wave  number,  3 

A  0  ~  o 

is  the  equilibrium  pinch  radius,  and  vA  is  the  Alfven  speed.  The  early 

experiments  were  indeed  observed  to  be  strongly  unstable,  in  general 

agreement  with  this  picture.  These  conclusions  have,  over  the  years, 

suppressed  interest  in  the  Z-pinch  as  a  quiescent  plasma  confinement 

geometry,  and  apparently  have  also  discouraged  theoretical  investigation  of 

the  stability  problem  from  a  more  sophisticated  point  of  view,  e.g.,  with 

the  inclusion  of  dissipative  and  kinetic  effects,  as  well  as  those  effects 

arising  from  time  dependence  of  the  equilibrium  state. 

However,  recent  experiments  with  deuterium-fiber-initiated  Z-pinches 
1  0 

at  NRL  and  Los  Alamos"  exhibit  apparent  stability  over  a  very  large  number 

of  Alfven  growth  times  a  /v  .  For  example,  in  the  NRL  experiment,  tne 

pinch  was  apparently  stable  over  the  entire  duration  of  the  current  rise, 

which  is  typically  130  ns  but  has  been  varied  down  to  110  ns  to  test 

sensitivity  to  this  parameter.  The  pinch  then  goes  sausage  unstable  after 

peak  current  is  reached.  Robson^  has  suggested  that  the  stability  is  a 

direct  consequence  of  rising  current  -  the  "I  hypothesis".  These 

observations  have  prompted  a  number  of  theoretical  efforts  to  re-examine 

stability,  with  the  emphasis  being  on  the  sausage  mode,  since  the  kink  mode 

is  not  seen  at  all  in  the  recent  experiments.  These  recent  and  ongoing 

efforts  include  linearized  one-dimensional  (radial)  treatments  which  assume 

a  mode  structure  exn(ikz)  a"d  model  fhe  plasma  cither  with  resistive  MHD.^ 

4 

or  MHD  with  the  full  Braginskii  transport  coefficients,  or  Chew- 

Manuscnpt  approved  August  14,  1989 
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Goldberger-Lov,  ’  or  Vlasov  ions.  In  addition,  there  have  been  two- 

g 

dimensional  nonlinear  fluid  simulations  that  include  resistivity  or  full 

4 

Braginskii  transport,  or  even  the  ongoing  ablation  of  the  deuterium 

9 

fiber.  Preliminary  results  from  some  of  this  work  has  suggested  that 

3  6 

resistivity  could  be  a  significant  stabilizing  factor.  ’  All  of  these 

treatments  are  essentially  numerical,  and  do  not  provide  analytic  scaling 

lavs,  nor  direct  insights  into  the  precise  nature  of  stabilizing 

mechanisms.  Also,  these  stability  calculations  have  in  most  cases  assumed 

1/3 

time  dependence  (t/t  )  for  the  current  I(t),  which  simplifies  and 
restricts  the  problem  by  leading  to  a  time-independent  pinch  radius.'^’1''’ 

In  this  work  we  investigate  the  resistive  MHD  sausage  linear  stability 
problem,  using  an  alternative  approach  which  has  the  virtues  of  simplicity 
and  generality,  and  directly  provides  analytic  scaling  laws  as  well  as 
insight  into  physical  mechanisms.  We  develop  a  model  that  admittedly 
includes  some  features  that  are  not  strictly  correct,  but  which  we  expect 
to  be  at  least  reasonable.  The  model  is  motivated  by  those  features 
appropriate  to  long-wavelength  modes,  but  one  hopes  that  it  applies 
qualitatively  for  kaQ  up  to  order  unity.  The  principal  assumptions  are  as 
follows:  (i)  The  plasma  temperature  T(r,z,t)  is  assumed  to  be  independent 

of  r,  in  both  the  equilibrium  and  the  perturbation,  (ii)  Self-similar 
equilibria  are  investigated  with  I(t)  «  (t/t  )a,  where  a  can  be  any  number. 

(iii)  No  heat  transport  along  z  is  included  in  the  investigations  to  date 
(but  this  could  easily  be  included,  and  would  be  destabilizing) - 

1?  13 

(iv)  Radial  pressure  balance,  i.e.,  the  Bennett  pinch  condition,  is 

assumed  in  both  the  equilibrium  and  the  perturbation,  (v)  The  current 

I(z,t)  is  assumed  to  be  unchanged  by  the  perturbation,  and  thus  independent 

2  2 

of  z.  (vi)  Some  terms  of  order  k  aQ  are  dropped  for  convenience.  (vii) 

The  fluid  motion  of  the  sausaging  is  assumed  to  be  self-similar.  However, 


2 


the  magnetic  field  B(r,z,t)  is  not  assumed  to  be  self-similar ,  but  rather 
is  calculated  from  resistive  MHD.  With  these  assumptions,  ve  are  able  to 
obtain  an  analytic  solution  for  the  grovth  of  sausage  perturbations. 

In  undertaking  this  calculation,  we  had  in  mind  two  particular 

resistive  effects  which,  coupled  to  increasing  I(t),  could  be  stabilizing. 

The  first  is  that,  according  to  the  Bennett  relation,  the  necks  of  the 

sausages  must  be  hotter,  and  therefore  less  resistive,  than  the  bulges; 

thus,  if  I  >  0,  the  skin  effect  is  enhanced  in  the  necks  and  (under  some 

circumstances)  the  current  in  the  necks  is  concentrated  in  a  smaller 

fraction  of  the  pinch  area  at  the  outside.  This  would  reduce,  at  the 

necks,  the  magnetic  pressure  in  the  interior  which  pushes  plasma  axially 

out  of  the  necks  and  thus  drives  the  instability.  The  second  effect  is 

that  (under  some  circumstances)  ohmic  energy  dissipation  is  enhanced  at  the 

necks,  where  the  plasma  is  compressed,  and  reduced  at  the  bulges.  This  can 

be  interpreted  mathematically  as  an  effective  decrease  in  the  plasma 

compressibility,  i.e.,  an  increase  in  the  adiabatic  index  from  its  actual 

value  5/3,  as  regards  the  perturbation.  We  have  calculated  these  and  other 

resistive  effects  fully  self-consistently ,  within  the  model  described 

above.  We  find  that  if  a  <  1,  the  growth  rate  is  reduced,  but  not  by  more 

than  about  a  factor  of  two,  from  the  well-known  value  from  ideal  MHD  with 

14 

an  assumed  time-independent  equilibrium.  In  fact,  there  are  other 
resistive  terms  which  are  destabilizing  (as  was  recognized  by  R.  J. 

Tayler^  many  years  ago),  and  also  some  geometric  terms  associated  with  the 
time  variation  of  the  equilibrium  radius  if  a  *  1/3.  If  a  >  1,  the  net 
result  is  that  the  growth  rate  is  larger  than  the  value  from  ideal  MHD  with 
a  time-dependent  equilibrium.  Our  model  also  can  be  interpreted  as 
representing  a  plasma  with  large  shear  viscosity.  Thus,  it  appears  that 
within  the  context  of  a  linearized  fluid  treatment  neither  resistivity  nor 


3 


shear  viscosity  can  account  for  the  long  period  of  stability  observed  in 
the  deuterium-fiber  experiments. 


II .  Self-Similar  Equilibrium 

The  "equilibrium"  state  of  the  pinch  is,  in  general,  time-dependent, 

since  I ( t )  will  have  a  (specified)  time  dependence  and  the  pinch  will  heat 

up  resistively  and  expand  or  contract  to  maintain  radial  pressure  balance. 

However,  by  definition,  the  equilibrium  is  z-independent .  We  assume  radial 

pressure  balance  is  maintained,  and  that  the  equilibrium  is  isothermal  and 

has  equal  electron  and  ion  temperatures.  (However,  all  the  conclusions 

also  follow  if  T  >>  T.  or  if  T  /T .  is  a  fixed  ratio.)  We  also  assume  that 
e  1  e  1 


the  conductivity 


a  £  T 


3/2 


(1) 
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in  accordance  with  the  Spitzer  formula  with  the  variation  of  the  Coulomb 
logarithm  neglected,  so  that  <jq  is  also  uniform  within  the  plasma.  The 
continuity,  radial  momentum,  and  magnetic  diffusion  equations  are  then 


3p, 


-o  1  3  n 

g—  +  -  -r—  rv  p  =  0, 
3t  r  9r  ro  o 


(2) 


3P  ,  B 
_ o  3 _ o 

9r  +  3r  8n 


o 

4  nr ' 


(3) 


3B 
_ c 

3t 


i_  v  B  .  L  i  L  rB 

3r  ro  o  4na  3r  r  3r  o' 
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(4) 


where  p  is  the  density,  P  the  pressure,  the  radial  velocity,  and 
subscript-zero  denotes  equilibrium.  There  is  also  an  energy  equation  which 
we  will  write  down  later  in  radially-integrated  form. 
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T„.  .  10,11,16.,. 
It  is  well  known  that 

solution,  defined  by  the  relations 


these  equations  admit  a  self-similar 


x 


ao(t)’ 


(5) 


V 

ro 


(t)x, 


(6) 


N  m . 

p0(r’t)  -  -fr~  p0(x)-  <7> 

Rao<t) 


Vr'l)  =  Bo(X)’  (8> 

0 


Jo<r,t)  *  7^77 100  ■  <9> 

Mo(t) 

where  a  (t)  is  the  Dinch  outer  radius,  N  is  the  number  of  ions  per  unit 
o  o  r 

length,  and  riK  is  the  ion  mass. 

Equation  (2)  is  automatically  satisfied  by  (5)  -  (7).  Equation  (3) 
implies  the  3ennett  relation, 


T 

o 


2’ 

4N  c 
o 


(10) 


where  the  temperature  T  is  in  ergs.  Beyond  that,  we  shall  not  need  to 
solve  (3)  for  the  pressure  P  .  Equation  (4)  reduces  to  a  Bessel  equation, 


I-  I  I-  x  B 

oX  X  dX  0 


4  rtcr  a  •  . 

O  O  In  n 
~  i  B0  = 


(11) 


to  which  the  solution  [normalized  to  unity  at  x  =  1  as  required  by  (8)]  is 


Bq(x)  = 


i-^qx) 

•  i \  ( q ) 

-i 

(  X 


,  0  <  x  <  1 

,  1  <  X 


(12a) 

(12b) 
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(13) 


where  is  a  modified  Bessel  function  and 


2 

q 


4ix<y  ( t  )a  ( t ) 
o  o 


2 

c 


till 

Kt)  * 


However,  we  have  assumed  B  (x),  and  therefore  q,  to  be  time-independent. 
This  can  be  satisfied  if 


I(t)  =  I0(t/to)a, 


in  which  case  (1)  and  (10)  imply 


ao(t)  «  (t/to)J“, 


and  thus  it  is  required  that 


aQ(t)  -  <t/tQ) 


(l-3a)/2 


so  that 


(14) 


(15a) 


(15b) 


2 

q 


a 

t 


4  rtcr  a 
o  o 


(16) 


is  indeed  time-independent. 

The  energy  equation,  including  both  ohmic  and  compressional  heating 
and  assuming  an  ideal  gas  with  adiabatic  index  5/3,  is 


3 


NT  t 


o  o 


4a 

—  NT 
a  oo 
o 


dr  2nr  J^. 

o 


(17) 


Using  the  Bennett  relation  (10)  to  eliminate  T  ,  Eqs.  (9),  (12)  and  (13) 
determine  JQ(r,t),  and  Eqs.  (14)  and  (15b)  to  specify  i/I  and  aQ/ao,  Eq . 
(17)  reduces  to  an  equation  that  determines  q  as  a  function  of  a: 


to 
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4a  = 


(18) 


I:2  (q> 

J  dx  x  I"(qx) 

0  0 

Equation  (18)  is  plotted  in  Fig.  1.  In  the  limiting  cases  of  small  and 
large  a,  the  Bessel  functions  can  be  expanded  and  ve  find 


q 


a  <  1 
a  »  1 


(19a) 

(19b) 


Also , 


in  the  case  a  =  1/3  which  corresponds 
q(o  =  1/3)  =  1.65 


to  a  fixed  radius  a  , 

o 

(19c) 


Recently,  Rosenau,  Nebel  and  Lewis^  have  shown  numerically  that  the 
self-similar  solutions  are  stable,  and  are  attractors  for  a  broad  class  of 
initial  conditions,  as  long  as  the  pinch  is  legislated  mathematically  to  be 
uniform  in  z. 

In  this  paper,  we  shall  apply  our  perturbation  model  to  this 
particular  equilibrium.  However,  the  perturbation  analysis  is  more 
general,  and  could  be  applied  to  other  assumed  equilibria.  Ve  have  in 
fact,  as  a  check  on  the  perturbation  model,  used  it  to  calculate  ideal  MHD 
growth  rates  for  various  assumed  time-independent  equilibria,  but  these 
calculations  will  not  be  presented  here. 
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Ill .  Linear  .zed  Perturbation  Analysis 

In  the  present  section,  ve  allow  all  quantities  to  depend  on  z  as  well 
as  r  and  t,  except  that  we  specify  that  I ( t )  is  z-independent  so  as  to 
prevent  charge  accumulation  anywhere  in  the  pinch.  Ue  shall  soon 
specialize  to  the  case  of  small  departures  from  a  z-independent 
equilibrium,  but  we  begin  by  writing  down  the  nonlinear  equations  for 
continuity,  momentum,  energy,  and  magnetic  field  diffusion.  Vith  long 
wavelength  perturbations  in  mind,  we  continue  to  assume  that  radial 
pressure  balance  is  satisfied,  and  that  T(r,z,t)  is  independent  of  r  but 
that  there  is  no  heat  conduction  a)ong  z.  The  basic  equations  then  are: 


i £ 
at 

*  5J  <ovz) 

1  3 

+  r  3r 

2 

2 

a 

B 

3r 

lP  8  nj 

4rir  ’ 

=  0, 


3v  *  /  n2. 

» tr  -  -  h  (p  *  y ' 


9B  a  f  3  f  c  a  1  a 
3t  +  3r  /r  J  +  3z  lz  J  3r  r  3r  r 


ci  a_  1  3B 
Art  3z  a  az' 


(20) 

(21) 

(22) 

(23) 


Now  we  make  the  key  simplification  that  allows  us  to  develop  a 
solvable  model:  we  require  the  fluid  motion  to  be  radially  self-similar, 
i.e.,  that  p(r,z,t)  be  of  the  form 

N(z,t)m. 

p(  r ,  z ,  t )  =  — 2 - -  P(x-Z)-  (24a) 

rta  (z,  t) 

where  N(z,t)  is  the  perturbed  number  of  ions  per  unit  length,  a(z,t)  is  the 
perturbed  outer  radius  of  the  pinch,  and  now 

x  5  r/a(z , t ) .  (25) 
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This  is  equivalent  to  specifying  that 


vr(r,z,t)  =  x  a(z , t ) ,  (24b) 

v  (r,z,t)  =  vz(z,t),  independent  of  r.  (24c) 

We  have  seen  that  self-similar  forms  vere  exact  solutions  for  the 
equilibrium,  but  Eqs.  (24)  are  not  in  general  exact  solutions  of  the  z- 
dependent  perturbation  equations;  they  are  imposed  externally  to  specify 
our  model.  Hovever,  ve  regard  them  as  a  qualitatively  reasonable 
representation,  particularly  in  the  long-wavelength  limit,  which  is  the 
focus  of  our  study,  for  the  class  of  sausage  perturbations  that  are  readily 
observed  experimentally  and  that  are  expected  to  be  most  disruptive  -  i.e., 
those  that  consist  of  a  more  or  less  coherent  pinching  of  the  entire  plasma 
cross-sec.tion  at  the  "necks",  and  expansion  of  the  entire  cross-section  at 
the  bulges.  Our  model  excludes  modes  which  are  strongly  localized  to  a 
limited  shell  in  r,  and  higher  radial  modes  which  may  have  parts  of  the 
plasma  cross-section  moving  inward  at  the  same  time  when  other  parts  are 
moving  outward.  In  some  cases,  particularly  within  ideal  MHD,  these  types 
of  modes  may  have  the  largest  linear  growth  rates,  but  one  does  not  expect 
them  to  have  the  same  type  of  disruptive  effect,  when  they  reach  large 
amplitude,  as  the  global  sausage  modes.  Since  we  are  excluding  certain 
classes  of  modes,  in  general,  ve  expect  our  model  to  yield  linear  growth 
rates  which  ore  so- jvhat  less  unstable  than  the  exact  results  of  resistive 
MHD,  buc  '  yield  qualitatively  accurate  growth  rates  for  any  global  modes 
that  are  prn'-rt.  Indeed,  we  have  tested  the  model  by  applying  it  with 
zero  resistivity  to  some  simple  time-independent  equilibria,  and  found  this 
to  be  the  case. 
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It  is  also  worth  noting  that  Eq .  (24c),  which  is  equivalent  in  itself 


to  the  self-similar  assumption,  becomes  exact  in  the  limit  of  long 
wavelength  and  large  shear  viscosity.  Thus,  the  present  model  becomes  a 
true  representation  of  the  physics  in  that  limit. 

The  inexactness  of  the  model  becomes  clear  when  Eq.  (22)  is  solved;  in 
general,  the  solution  will  be  of  a  value  of  v^  that  varies  with  r, 
contradicting  Eq.  (24c).  Thus,  to  make  our  model  self-consistent,  we  set 


v^  equal  to  the  average  value  given  by  (22),  i.e.,  we  replace  (22)  with 


Now  we  linearize  our  basic  set  of  Eqs.  (20),  (21),  (23),  (28),  and 
(30),  by  defining 


aQ(t)  - 

(31a) 

<J(  z .  t )  ,  ,  * 

»„<'>  ’ 

(31b) 

N(zr 12  -  1  +  N 

N0(t)  ~ 

(31c) 

I(t)  =  IQ(t), 

(31d) 

B(  r ,  z ,  t )  =  ~  B0  (~  r,t)  +  ^  B  ( x ,  z ,  t ) , 

(32a) 

J(r,z,  t)  =[^j  JQ  r»t  )  +  J(x>z>t)- 

^  ~  fta 

(32b) 

Equations  (32)  define  the  first-order  quantities  B  and  J  in  terms  of  the 
difference  between  the  exact  values  B(r,z,t),  J(r,z,t)  and  self-similar 
forms  that  carry  the  same  current.  This  somewhat  unusual  approach  is  very 
convenient  for  handling  the  boundary  conditions  at  r  =  a,  where  both 
J(r,z,t)  and  <j(r,z,t)  jump  discontinuously  to  zero.  With  the  definitions 
(32)  and  using  x  rather  than  r  as  our  independent  variable,  all 
discontinuities  occur  exactly  at  x  =  1,  it  is  not  necessary  to  add  any 
surface  terms  to  the  first-order  equations,  and  the  boundary  conditions  are 

B(0,t)  =  0,  (33a) 

B(x, t)  =  0,  for  x  >  1,  (33b) 

J(x,t)  =  0,  for  x  >  1.  (33c) 

Note  that  J  is,  in  general,  discontinuous  at  x  =  1,  but  if  a  is  finite  B  is 
continuous . 
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Ue  treat  all  quantities  with  carets  as  first-order  small  quantities 
with  z-dependence  exp(ikz),  neglect  second  order,  and  also  define 


v  =  ikv  . 
z 


(34) 


Then  Eq.  (2.0)  reduces  to 


dN 

dt 


=  -v . 


(35) 


Equations  (28),  (31d),  and  (35)  give 


(36) 


3/2 


and  continuing  to  use  a  *  T  ,  we  have 


da  3  * 
dt  =  2 


(37) 


It  is  useful  to  explicitly  expand  Eqs.  (32)  to  first  order,  using  (8) 
and  (9) 


B( 


r,t)  =  ~  [d  -  a)  Bo(x)  +  B(x,t)], 


J(r , t)  =  [(1  -  2a)  J0(x)  +  J(x*t)]- 

na 

o 


We  integrate  (21),  using  (38),  (12b),  and  (33b),  to  obtain 


(38a) 


(38b) 


p  +  _  I*  dr'B2(r') 

811  ”  J  4 nr' 


1  , 


„  2  2 
nc  a 

o 


[(1  -  2a)  {j  +  J  ~~  B2(x')  +  2j  Bq(x')B(x')] 

x  x  X  x 

(39) 
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Ve  now  use  this  in  the  equation  of  motion  (26).  Reversing  the  order  of 
integrations  and  explicitly  using  (12)  for  the  equilibrium  field,  ve  obtain 


2  2 

dv  2k  I 

" 

f  i 

.  dxxl^(qx)!  „  dxxI1(qx)B(x)' 

dt  -  M  2 

N  m .  c 

O  1 

-a 

[  2  +  J 

^  C 

,  J  ‘  1  ^  J 

(40) 


After  some  algebra  but  no  further  approximations,  the  magnetic  field 
Eq.  (23)  can  be  put  in  the  form 


9_ 

at 

v 


a 

t 


4na  a 


3_ 
2  3x 


1  3_ 
x  3x 


•x  + 


,2  2 

k  c 


4rta 


oj 


a  ‘  a  ,  ‘  0  \ 

-  a  -  v  -  -  (a  +  3a) 


k"c" 

4rta  3 
o 


I2(qx) 


(41) 


Finally,  the  ohmic  heating  term  in  (30)  can  be  calculated  in 
linearized  form  and  the  energy  equation  put  in  the  form 


da 

dt 


5  * 
4  V 


l-3« 
2 1 


A  - 


4tta  a^ 
o  o 


'  2  . 


-  2a)  +  4q^ 


dxxl^ (qx)B(x) 


Ix(q) 


vhere  we  have  explicitly  used  the  time  dependence  (15b)  of  aQ(t). 

Equations  (37),  (40),  (41),  ana  (42)  with  the  boundary  conditions 

(33a, b)  constitute  a  complete  set  of  homogeneous  equations  for  the  time 

evolution  of  <x(t),  v(t),  a(t),  and  B(x,t).  Ve  clean  up  the  equations  a  bit 
2  n 

by  dropping  the  k^c4"  terms  on  each  side  of  (41);  these  terms  correspond  to 
magnetic  diffusion  along  z,  and  are  stabilizing  but  small  in  the  long 
wavelength  limit.  We  also  use  (13)  and  (14)  to  eliminate  in  terms  of  q, 
and  define  a  characteristic  Alf^en  velocity  v,  by 


2 

B  (a  ) 
o  o 


(  21 
4 n  N  m. /rta 
V  O  1  o) 


N  m .  c 
o  1 


2 ' 


(43) 
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Then  Eqs.  (40)  -  (42)  become 


dv  .2  2 
j-  =  k  v 
dt  A 


[1  +  2In(q)I  a  +  2 


dxxI^(qx)B(x, t ) 


IjCq) 


5  l-3a  *  a  a  4a  T  dxxI1(qx)B(x, t) 

4  V  2t  a  '  t  2 1  t  J  I. (q)  ’ 

0  1 

-  f  -  b  *  h  x)B(x>t)  =  [“  V  "  ?  (a  +  2a)]  i^T 


,(46) 


where 


In(q)  - 


dxxl- (qx) 
n  * 


^(q) 


/  t  ~1  ^  \ 


The  quantity  1^,  which  can  be  regarded  as  a  function  of  a,  is  shown  in 
Fig.  1.  Using  (18)  and  doing  some  integrations  by  parts,  we  can  write 


yq) 


■11  -  qyq)  ~  4a' 


(47b) 


in  which  form  it  is  clear  that 


Ill(0)  =  4’  if  a  << 


I^^(®)  =  if  a  >>  1. 


(48a) 


(48b) 


1^  is  a  monotonic  decreasing  function  of  a  (i.e.,  of  q)  which  is  well 
approximated  by 


hlM  =  Ml^o*  (49) 

The  expression  (49)  is  always  smaller  than  the  exact  value  of  1^,  with  the 
inaccuracy  peaking  at  -  12%  over  a  wide  range  of  a  from  -  0.5  to  ~  3. 


Before  proceeding  with  the  solution  of  Eqs.  (37),  (44)  -  (46),  and 
(33),  ve  shall  pause  to  discuss  the  general  nature  of  these  equations  and 
the  origin  and  significance  of  their  terms.  If  instability  is  present,  in 
general,  all  time  derivatives  will  be  positive  at  times  of  primary  interest 
(instability  is  non-oscillatory ) ,  as  vill  be  discussed.  At  a  sausage  neck, 
vhere  a  <  0,  ve  vill  find  v  >  0  (plasma  floving  out  axially),  and  thus,  N  < 
0,  and  thus  a  >  0  (the  neck  must  be  hotter,  according  to  the  Bennett 
relation,  and  so  has  higher  conductivity).  The  basic  instability,  vhich  is 
familiar  from  ideal  MHD  analysis  of  a  time-invariant  equilibrium,  is  given 
by  the  first  terms  on  the  right-hand  side  (RHS)  of  (44)  and  (45);  these 
yield  a  growth  rate  on  the  scale  of  kv^,  vith  some  dependence  on  the 
equilibrium  current  profile  through  1^. .  The  last  term  of  (44)  represents 
the  effect  of  non-self-similar  evolution  of  B:  typically  J(r,z,t)  vill  be 
more  peaked  at  the  outside  of  the  pinch  at  a  neck,  i.e.,  B(x)  for  0  <  x  <  1 
vill  be  negative  vhere  a  <  0,  vhich  veakens  the  pressure  pushing  plasma  out 
of  the  neck.  Thus,  this  term  is  stabilizing.  In  Eq.  (45),  the  term  (1- 
3a)a/2t  arises  simply  from  the  fact  that  perturbations  to  a(z,t)  must  be 
scaled  to  aQ,  vhich  is  itself  a  time-varying  quantity;  this  geometric 
effect  does  have  a  significant  effect  on  instability  grovth  rates,  as  we 
shall  see.  The  last  three  terms  of  (45)  arise  from  the  differences  in 
ohmic  heating  between  sausage  necks  and  bulges:  if  there  is  extra  heating 
at  the  necks,  this  is  stabilizing  by,  in  effect,  increasing  the  adiabatic 
index  (stiffness)  of  the  plasma  and  opposing  further  shrinkage  of  the  neck. 
The  a/t  term  in  (45)  represents  increase  of  ohmic  heating  at  a  neck  due  to 
self-similar  compression  of  the  current  density;  the  a/2t  term  represents 
decrease  of  ohmic  heating  due  to  increased  conductivity;  and  the  last  term 
represents  increased  ohmic  heating  due  to  current  concentration  at  the 
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outside  resulting  from  departures  from  self-similari ty .  It  remains  to  be 
determined  under  what  circumstances  the  stabilizing  terms  win  out. 

Equation  (46)  shows  the  ways  in  which  B(r)  departs  from  self¬ 
similarity.  The  first  term  on  the  RHS  represents  decrease  of  B  at  a  neck 
due  to  axial  convection  of  B  along  with  the  outflowing  plasma.  The  a(  a  + 
2a)/t  term  represents  the  first-order  contribution  to  the  skin  effect  in  a 
situation  where  I  >  0.  As  in  the  discussion  of  Eq.  (45),  the  strength  of 

the  skin  effect  scales  as  the  ratio  of  the  magnetic  diffusion  time 
2  2 

4 rtcr  a  /c  to  the  current  rise  time  I/I;  hence,  the  first-order  term  scales 
oo 

as  a(  cr  +  2a)/t.  Clearly,  these  terms,  which  arise  from  the  combined  effect 
of  resistivity  and  rising  current,  have  the  potential  to  be  stabilizing, 
but  must  be  evaluated  to  determine  their  effect. 


IV.  Solution  of  the  Perturbation  Equations 

Equations  (37),  (44)-(46),  and  (33a, b)  could  be  solved  numerically, 
e.g.,  by  using  a  Fourier-Bessel  expansion  of  (46).  However,  as  our 
objective  here  is  to  obtain  analytic  scaling  laws  in  simple  form,  we  choose 
to  perform  a  simpler  approximate  analysis.  We  note  that  B(x,t)  appears  in 
Eqs.  (44)  and  (45)  only  in  the  integral  form 


1 


b(t)  *  I 


0 


dxxI^(qx)B(x, t) 


ix(q) 


Equation  (46)  for  B  has  the  particular  solution 


(50) 


B1(x, t)  =  F(t) 


I1(qx) 

f^TqT’ 


(51) 


F(t)  s  J  dt'  v  +  “(cr  +  2  a) 


(52) 
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However,  (51)  fails  to  satisfy  the  boundary  condition  (33b),  and  thus,  we 
must  add  in  a  homogeneous  solution  B^ix.t)  with  the  correct  boundary 
condition: 

B(x,t)  =  B1(x,t)  +  fi2(x,t),  (53) 


a 

t 


a 

q2t 


h  ?  h  x)b2  -  °- 


B2(0,t)  «  0, 
B2(l,t)  =  — F ( t ) . 


(54a) 

(54b) 

(54c) 


Ue  shall  argue  that  over  the  time  scale  of  primary  interest  to  us,  it  is 
reasonable  qualitatively  to  neglect  B2- 

In  general,  we  cannot  produce  a  closed-form  solution  for  B2(x,t).  But 
for  the  special  case  where 

F(t)  -  tn,  (55) 

the  solution  to  Eq.  (54)  is 

I.  ( /1+q/a  qx) 

B„(x,t)  =  -F(t)  -  -  - - .  (56) 

I1(/l+n/a  q) 


In  this  case,  the  integral  b(t)  of  Eq.  (50)  can  be  written  as  a  sum  of 
contributions  from  B^(x,t)  and  B2(x,t): 


b(t)  =  bx(t)  +  b2(t), 


bj(t)  =  F( t) 


1 

dxx 

0 


I2(qx) 

lj(q) 


(57) 


(58a) 


.  I. (qx)I. (/1+q/o  qx) 

b-(  t)  =  -F(  t )  dxx  — -  ;  ■  - .  (58b) 

0  I1(q)I1(/l+n/«  q) 
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If 

^  q2  »  1  ,  (59a) 

and  also 

^  q2  »  q2,  (59b) 

then  Bn(x,t)  is  locali2ed  near  x  =  1  and  as  a  result 

b2  >>  br  (60) 

Conditions  (59)  are  both  satisfied  if  and  only  if 


Hq *" / a  > >  1 .  (61) 

2  1/2 

In  fact,  b^/b^  is  of  order  (a/qq  )  .  The  interpretation  of  these 

equations  is  that  if  the  perturbation  grows  rapidly  enough  compared  to  the 
equilibrium  current-rise  time  scale,  then  B2(x,t)  is  essentially  a  surface 
current  perturbation,  which  does  not  have  time  to  diffuse  substantially 
into  the  pinch;  this  current  contributes  little  to  the  integrated  quantity 
b(t)  that  occurs  in  (44)  and  (45). 

If  the  time  dependence  of  F(t)  is  not  a  power  law,  we  cannot  solve  for 
B?(x,t)  in  closed  form,  but  the  same  physical  argument  applies  with  h 
defined  by 


n  ,  SLS 11. 

n  F(t) 


(62) 


the  condition  for  neglecting  is 

(63) 


For  the  limiting  cases  of  small  and  large  a,  we  see  from  Eq.  (19)  that  this 
reduces  to 
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(64a) 


yt  »  |  ,  if  a  «  1, 
or 

yt  »  1/ 4a,  if  a  >>  1.  (64b) 

In  either  case,  the  criterion  is  well  satisfied  at  a  time  when  no  more  than 
one  e-fold  of  unstable  growth  has  occurred.  We  are  generally  interested  in 
the  longer  time  scale  over  which  an  instability  e-folds  several  times,  so 
it  is  reasonable  to  approximate  B(x,t)  by  the  particular  solution  F^(x,t) 
whose  x-dependence  is  known.  Ve  note  from  (51),  (56)  and  (58)  that  this 
approximation  increases  |b(t)|,  and  this  is  true  quite  generally  if  I  >  0. 
The  quantity  b(t),  i.e.,  the  effect  of  the  deviation  of  B(r,t)  from  self¬ 
similarity,  is  unambiguously  stabilizing,  as  will  be  discussed  below. 

Thus,  the  approximation  of  B(x,t)  by  B^(x,t)  reduces  the  growth  rates  which 
we  shall  calculate. 

Referring  to  Eqs.  (50)  -  (52),  (57)  and  (58),  ve  see  that  the  partial 
differential  equation  (46)  can  now  be  replaced  by  an  ordinary  differential 
equation  for  b(t),  so  that  (37),  (44)  -  (46),  and  (33a, b)  reduce  to  a  set 
of  four  homogeneous  first-order  linear  ordinary  differential  equations, 
which  can  be  written  in  the  matrix  form 

v  1  ro  (U2I..,  )k2v2  -  2k2v2  0  1  fv 

11  A  A 

a  5/4  3(l-a)/2t  4a/t  l/2t  a 

b  2I^a/t  0  I^a/t  b 

a  J  [-3/2  0  0  0  a 

The  eigenvalues  of  this  matrix  give  the  instantaneous  growth  rate  of 
perturbations;  this  growth  rate  will  generally  be  time-varying,  since  the 
equilibrium  varies  with  time,  i.e.,  growth  wi_x  not  be  purely  exponential. 


19 


It  is  convenient  to  define  T  as  the  eigenvalue  scaled  to  kv^(t).  We  then 
find  the  following  quartic  equation  to  determine  the  four  eigenvalues: 

0  =  4(kvAt)2T4  +  6(l-a)kvAt  T3  -[(5+2In)(kvAt)2  +  32a2Injr2 

+  [_3+(6-4a)In  +  32<xIn2jkvAt  T  +  6(l+a+8oI11)aI11 .  (66) 

The  ideal  MHD  limit  corresponds  to  kvAt  -*  (The  temperature 

2 

increases  steadily  as  T  11  I  ,  so  resistivity  becomes  negligible  at  late 
times.)  It  is  also  possible  to  express  kvAt  as 

kv  t  =  (ka  )S,  (67) 

A  10 

q 


4naQ(t)a2(t)  vA(t) 

7 


(68) 


18  3  8 

is  the  Lundquist  number,  which  has  been  used  in  the  literature  ’  as  an 

indication  of  the  strength  of  resistive  effects.  It  is  apparent  from  Eqs. 

(65)  -  (67)  that  the  growth  rate  and  other  characteristics  of  the 

instability,  for  a  given  value  of  a,  can  depend  on  S  and  kaQ  only  as  the 

combination  (kaQ)S.  S  is  a  monotonically  increasing  function  of  time,  and 

large  S  corresponds  to  the  ideal  MHD  limit.  In  this  limit,  (66)  gives 


i.e.,  growth  is  on  the  Alfven  time  scale,  as  expected.  Ue  may  regard  the 
four  eigenvalues  T^,  and  F^  as  functions  of  kvAt.  The  eigenvalue 

that  appears  in  Eq.  (69),  which  we  shall  denote  F^(kvAt),  is  always  the 
largest  and  thus  of  primary  interest.  Ue  have  solved  Eq.  (66)  for  the 
value  of  T  as  a  function  of  kvAt  or  equivalently  of  (kaQ)S.  It  is 
convenient  to  plot  as  a  function  of  (kv  t)-^1,  which  is  an  indication  of 
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the  strength  of  non-ideal  effects.  The  results  are  shown  in  Fig.  2.  It  is 
apparent  that  for  a  <<  1,  the  value  of  is  reduced  by  non-ideal  effects, 
but  only  to  =  1/2  even  for  a  -»  0  in  the  limit  kv^t  •*  0.  For  small 
finite  values  of  a,  e.g.  the  case  a  =  1/3  shown  in  Fig.  2,  F^  decreases  as 
a  function  of  (kv  t)  \  but  goes  complex  for  times  earlier  than  a  finite 
value  kv.t*.  Thus,  the  instability  is  predicted  to  exhibit  oscillatory 
growth  for  t  <  t*.  However,  t  <  t*  corresponds  to  less  than  one  e-fold  of 
growth,  and  our  neglect  of  B?(x,t)  [Eqs.  (53)  -  (58)  J  becomes  questionable 
at  these  early  times,  so  ve  do  not  plot  the  complex  values  of  F^ .  In  no 
case  do  ve  find  stability,  or  even  reduction  of  the  ideal  MHD  growth  rate 
by  more  than  aoout  a  factor  of  two,  for  time  scales  of  interest. 

For  a  >  1,  ve  find  to  be  a  rapidly  increasing  function  of 

(kv^t)  1,  i.e.,  the  growth  rate  is  much  larger  than  the  ideal  MHD  growth 

rate  for  a  static  equilibrium.  Tracing  back,  through  the  mathematics,  we 

can  see  that  this  rapid  growth  is  primarily  due  to  the  2-2  matrix  element 

of  Eq .  (65),  i.e.,  to  the  term  (l-3a)a/2t  on  the  RHS  of  Eq.  (45)  for 

da/dt.  This  is  not  a  resistive  effect,  but  simply  a  geometric  one  arising 

from  the  increasing  I(t),  which  is  not  oaianced  by  enougn  onmic  heating  to 

maintain  pressure  balance  at  constant  pinch  radius  aQ.  If  a  >  1/3,  the 

equilibrium  radius  at  t)  of  the  pinch  steadily  decreases  as  t^i_  a  '  “ ,  but 

any  perturbations  to  the  radius  do  not  collapse  as  quickly.  This  effect 

tends  to  increase  a(t)  =  (a  -  a  )/a  ,  which  is  the  quantity  that  drives  the 

o  o 

growth  of  the  perturbations.  Additionally,  because  of  this  effect 

2 

perturbations  to  aa  are  dominated  by  the  2a  perturbation  rather  than  the  a 
perturbation:  hence  the  skin  effect  is  weaker  at  a  sausage  neck,  which  is 
destabilizing,  as  discussed  in  Sec.  III. 
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Concluding  Remarks 


Ue  have  used  a  self-similar  model  to  calculate  sausage  mode  growth 

rates  within  the  framework  of  resistive  MHD.  Although  our  model  is  a 

simplification  of  the  perturbation  dynamics,  it  does  treat  quite  a  general 

class  of  time-dependent  exact  equilibria,  it  is  fully  self-consistent  in 

its  treatment  of  equilibrium  and  perturbations,  and  it  is  a  fully  resistive 

treatment  from  the  outset  rather  than  being  a  perturbative  treatment  of 

departures  from  ideal  MHD.  Furthermore,  we  expect  the  self-similar 

assumption  to  be  at  least  qualitatively  reasonable  in  calculating  growth 

rates  for  global  modes,  and  in  any  event  to  yield  growth  rates  that  are 

smaller  than  exact  solutions  of  resistive  MHD.  Thus,  it  is  disappointing 

that  our  results  fail  to  explain  the  long-lived  apparent  stability  of  the 

1  2 

recent  deuterium-fiber  Z-pinch  experiments.  ’  Ve  do  find  reductions  in 

1/3 

the  instability  growth  rate  for  the  case  I  «  t  where  the  equilibrium 
pinch  radius  is  time-independent,  and  for  other  cases  I  «  ta  with  a  small, 
but  the  reductions  are  by  no  more  than  a  factor  of  two,  which  is 
insufficient  to  explain  the  experiments. 

Assuming  the  experimental  observations  to  oe  a  valid  indication  of 
stability  over  a  significant  parameter  regime,  we  may  well  ask  what  types 
of  physical  effects  could  account  for  this  stability.  Several  types  of 
effects  come  to  mind,  which  have  been  excluded  by  the  assumptions  of  the 
present  model,  and  also  by  assumptions  made  in  other  recent  calculations. 
First,  we  (and  others)  have  assumed  the  current  I(z,t)  to  be  z-independent 
and  equal  to  the  equilibrium  current.  If,  in  fact,  I(z,t)  is  reduced  at 
the  sausage  necks  -  which  would  amount  to  some  capacitive  loading  with 
charge  accumulating  on  the  shoulders  of  the  sausages  -  that  would  have  a 
stabilizing  effect  on  the  mode.  Secondly,  we  (and  others)  have  assumed  the 
temperature  within  a  given  slice  (i.e.,  at  a  given  value  of  z)  to  be 
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independent  of  r.  This  is  in  agreement  with  the  numerical  equilibrium 

1  o 

calculations  of  Lindemuth,  et  al  provided  there  is  no  remaining  residual 
solid  fiber.  However,  these  authors  argue  that  many  tens  of  nanoseconds 
pass  before  the  fiber  is  fully  ablated,  and  find  that  the  presence  of  the 
fiber  at  r  =  0,  which  cools  the  plasma  there,  results  in  a  ramped  profile 
for  T(r).  It  is  possible  that  such  a  profile  is  more  stable  to  sausage 
perturbations.  Ve  have  not  calculated  either  of  these  effects  and  merely 
suggest  them  as  possible  avenues  for  future  investigation. 

More  generally,  it  should  be  noted  that  diagnostics  used  in  the 
present  generation  of  fiber  z-pinch  experiments  would  detect  sausage 
instability  only  when  it  has  reached  large  amplitude.  Thus,  it  is  possible 
that  nonlinear  effects  excluded  from  the  present  linearized  treatment,  and 
possibly  coupled  to  the  effects  of  increasing  current,  resistivity,  and 
other  dissipative  mechanisms  may  play  a  role  in  explaining  the  stability 
observed  in  the  experiments. 
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a  =  ^ 


a  =  1/3 


Fig.  2.  The  instability  growth  rate  T  [normalized  to  kv^,  where  the  Alfven 
speed  vA  is  defined  in  Eq.  (43)],  for  four  different  values  of  the 
parameter  a,  which  specifies  the  time  dependence  of  the  current, 
that  the  value  of  T  for  a  =  1/3  becomes  complex  for  (kv  t)-^  > 
(kvAt*)_1  =  1.8. 
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